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Summary: We review and compare recent work on the properties of 
fluctuating interfaces between nematic and isotropic liquid-crystalline 
phases. Molecular dynamics and Monte Carlo simulations have been 
carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, 
and the fluctuation spectrum of interface positions (the capillary wave 
spectrum) has been analyzed. In addition, the capillary wave spectrum 
has been calculated analytically within the Landau-de Gennes theory. 
The theory predicts that the interfacial fluctuations can be described in 
terms of a wave vector dependent interfacial tension, which is anisotropic 
at small wavelengths (stiff director regime) and becomes isotropic at large 
wavelengths (flexible director regime). After determining the elastic con- 
stants in the nematic phase, theory and simulation can be compared 
quantitatively. We obtain good agreement for the stiff director regime. 
The crossover to the flexible director regime is expected at wavelengths of 
the order of several thousand particle diameters, which was not accessible 
to our simulations. 
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Introduction 

Elongated particles form liquid crystalline structures at high densities. For example, they 
often exhibit a nematic phase, where the particles align along one preferred direction (the 
director), while their positions are disordered, like in a fluid. For symmetry reasons, the 
phase transition between the nematic phase (N) and the fully disordered isotropic phase 
(I) must be first order. Hence the two phases coexist in a certain parameter regime, and 
are separated by a nematic-isotropic (Nl)-interface. The properties of that interface are 
intriguing for several reasons: Since it separates two fluid phases, the interfacial tension 



of a planar interface of given orientation N should be isotropic - the free energy cost of 
increasing the intcrfacial area should not depend on the direction in which the area has 
been extended. On the other hand, the interface interacts with the adjacent, anisotropic, 
nematic fluid - for example, it orients the director ("surface anchoring") - and one would 
expect this to have an influence on the interface. Hence, it should also exhibit anisotropic 
features. 

The role of the interfacial tension for the intcrfacial structure can be assessed by studying 
the fluctuations of the interface positions, the "capillary waves" . For interfaces between 
two simple fluids, the theory of capillary waves is quite simple: Let us consider a planar 
interface oriented in the z-direction with the projected area A, neglect overhangs, and 
parametrize the local interface position by a single- valued function h{x,y). Fluctuations 
enlarge the interfacial area and thus cost the free energyP"^! 



with the interfacial tension 7. Based on this functional, the thermal height fluctuations 
(|/i(q)P) at given wave vector q can be calculated in a straightforward manner. One 
obtains the capillary wave spectrum 



At a NI interface, the situation is more complicated. The capillary wave fluctuations 
are then determined by at least three factors, namely, (i) the interfacial tension, (ii) the 
surface anchoring, and (iii) the elasticity of the nematic bulk, i.e., the elastic response of 
the nematic fluid to local director variations. 

In the present paper, we shall review recent work that has shed light on the interplay of 
these three factors. The paper is organized as follows. In the next section, we discuss 
large-scale simulations of Nl interfaces in model liquid crystals, which have allowed to 
study in detail the capillary wave spectrum in these systems. In section three, a Landau- 
de Gennes theory of capillary waves is presented. The predictions of this theory are 
compared with the simulations in section four. Finally, we summarize and conclude with 
a brief outlook on nonequilibrium interfaces. 

Computer Simulations 

Computer simulations of equilibrium NI interfaces were carried out for two standard 
model liquid crystals: Systems of ellipsoids'^'^' and systems of spherocyhnders.'^"^' In 
two of these studies, '^''^l the system sizes were large enough that capillary waves could be 
investigated in detail. These shall be compared with each other. 

The model ellipsoids interact with each other via a soft, repulsive. Weeks- Chandler- 
Andersen-type[^°l potential 



J^{h} = 7 Jdx dy ^1 + + (dyh)^ ^7^ + 1 Jdx dy [{d^hf + [dyhf] (1) 
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where as is the diameter of the eUipsoids, Uj and \ij are the orientations of eUipsoids i 
and j, iij is the unit vector connecting the centers of the two eUipsoids, and the shape 
function^^^l 
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approximates the contact distance between hard eUipsoids of elongation k, with x = 
(k^ — + 1). The temperature was chosen k^T = e. 

The spherocyhnders are taken to be hard rods with semi-spherical caps, i.e., lines of length 
L, which may not come closer to each other than a distance D (the diameter of the rods). 
The elongation of the particles, k or L/D, respectively, was chosen 15 in both cases. The 
simulations were carried out in the A^V^T-ensemble with roughly ^ 100.000 particles in 
an elongated simulation box with side length ratios L^. : : = 1 : 1 : 2 and periodic 
boundary conditions. The simulation method was Molecular dynamics in the case of 
the ellipsoids,'^] and Monte Carlo in the case of the hard rods.'^^ By choosing a mean 
density between the densities of the nematic and isotropic phase at coexistence, phase 
separation into a nematic and an isotropic slab was enforced. The equilibrated systems 
thus contained two planar NI interfaces with orientation normals N parallel to the long 
axis of the simulation box. Both for ellipsoids and hard rods, the director in the nematic 
phase spontaneously aligned in the direction parallel to the interface (planar anchoring). 
Normal anchoring could be studied as weU, but had to be enforced, e.g., by choosing 
special boundary conditions, t^' or by preparing an initial configuration that contains a 
nematic slab with normal orientation. Here, we will almost exclusively discuss planar 
anchoring. To analyze the capillary wave spectrum, the simulation box was split into 
columns (blocks) B x B x Lz, and the local interface positions were determined separately 
in each block. This gave two interface topographies h{x, y) for every configuration, which 
could then be Fourier transformed to obtain the capillary wave spectrum, |/i(q)p (cf. 
Appendix). More details on the simulation and the data analysis can be found in Refs. 6 
and 9. 

Table [1] summarizes the main properties of our model systems at coexistence. Most of 
the data are compiled from earlier publications, l^'^'*^'^' but the table also shows new, 
previously unpublished results. In particular, it gives the Frank elastic constants'^' 
Ki, K2, of the ellipsoid system in the nematic phase, which we have evaluated in order 
to use them for a quantitative comparison between theory and simulation (see below). To 
this end, separate simulations of a homogeneous system at the density of the nematic phase 
were conducted, and the order tensor fluctuations were analyzed, following a procedure 
described in Refs. 13 and 14. 

The capillary wave spectra for the two model systems are shown in Fig. [TJ The capillary 
waves are clearly anisotropic: In the direction of the director, they are suppressed by a 
factor 0.3-0.5, compared to those in the perpendicular direction. Nevertheless, they still 
roughly follow the proportionality law 1 / {\h{c\)\'^) ~ predicted by Eq. ([2]). Deviations 
are observed in the ellipsoid system for wavelengths comparable to the particle length, 
27r/g ~ ?)K,as- In the hard rod system, the proportionality persists over the whole accessi- 
ble wave vector range. Hence it seems that Eq. ([2]) would provide a reasonable description 
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Table 1: Properties of our model systems at NI coexistence. Numbers in brackets indicate 
errors in the last digit. 





Ellipsoids 


Hard rods 


densities 


Pn = 0.0181(l)/crf 
Pi = 0.0149(l)/a3 


Pn = 0.027/L'3 
Pi = 0.023/D=' 


reduced densities 

* / 

P = PI pep 


P*N = 0.192 
p*j = 0.158 
Pep = V2/15 


Pn = 0.22 
p*j = 0.193 
Pep = 2/(V2 + 15V3) 


order parameter 


= 0.74151 


Sn ^ 0.7 


Frank 
elastic constants 
(nematic phase) 


= 0.78(2) 
K2 = 0.32(2) at p = 0.018 
Ks = 3.61(5) 




interfacial tension 
planar anchoring^ 

normal anchoring 


from pressure tensor'^'^' 


from histogram method'^' 


77JV = Ommk^T/ai^'i 
^iN = 0.009(3)A;^T/a2l6l 


^jj^ = 0.0064A;^T/D"^«I 


= 0.015(4)A;^T/or,^l5l 




interfacial width-f 


w ~ 7. Sag 


w ^ 10.6D 



^ The two values for the eUipsoids were obtained in simulations with different system sizes. They agree 
within the error. 

t The interfacial width has been determined by a fit of the local order parameter profile to a tanh profile, 
tanh((z — Z[))/w. Due to the effect of capillary waves, it depends on the lateral system size. The values 
given here correspond to local profiles evaluated in blocks of lateral size B — 2Kas (ellipsoids) and 
B^2L (rods). 

of our data, provided the interfacial tension 7 were allowed to depend on the direction of 
the wave vector, q. On the other hand, we have pointed out in the introduction that the 
macroscopic surface tension between two fluid phases must be isotropic. To resolve this 
apparent contradiction, a Landau-de Gennes type theory of capillary wave fluctuations 
at NI interfaces has been developed. '^^1 It shall be presented next. 

Theory 

The Landau-de Gennes theory is based on a free energy expansion in powers of a sym- 
metric and traceless (3x3) order tensor field Q(r)[^' 

F = JS {^Tr(Q2) + |Tr(Q3) + ^Tr(Q2)2 + h^,QJk^^QJk + Y^iQrAQkj}, (6) 

with B < close to the NI transition. If biaxiality can be neglected (both in our ellipsoid 
and hard rod system, the maximum value of the biaxiality near the NI interface was 
less than 0.06), the order tensor can be written as'^''! Qij{r) = |S'(r)(3?7,j(r)?T,j(r) — 6ij), 
where ^(r) is the local scalar order parameter, and n a unit vector characterizing the 
local director. For further simplification, we introduce "natural" units So = —2B/9C, 
^0 = 2y^(Li + L2/6)/3C ^o""*^, and eo = (3(7/16)5*0^0 for the order parameter, the length, 
and the energy, and rescale all quantities by these units. The free energy functional 
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Figure 1: Inverse mean-squared Fourier components of the interface position h vs. wave 
vector q squared for q parallel and perpendicular to the director ng. Left: Ellipsoids 
(similar to Ref. 6; data taken from Ref. 15). Right: Hard rods (similar to Ref. 9, different 
data set). The straight dashed lines indicate the prediction of Eq. ([2]) with the interfacial 
tension taken from Tabled], 7 = O.OOQ/c^T/ct^ (ellipsoids) and 7 = 0.0064A;^T/Z}^ (hard 
rods). Remarkably, it seems to fit the lower branch (q _L Hq) rather well in both systems. 

then takes the formt""^^! 

F = ^jd\{f + gy + g2} with f = S\{S-lf + t), (7) 

= (V-5)2 + a(n-V5)2 + 4a5((V-n)(n-V5) + ^(nx Vxn)(VS)), 
g^ = 5^(^(3 + 2a)(Vn)2 + (3 - a)(n ■ V x n)^ + (3 + 2a)(n x V x n)^' 

with the two dimensionless parameters, 

t=\A^-l, and a = l ^) (8) 
4 eo 2 (Li + L2/6) 

The parameter t gives the distance from NI coexistence and becomes zero at coexistence. 
The parameter a characterizes the response of the system to order parameter and director 
variations. In Eq. ([7]), the first term, f{S), is the free energy density of a homogeneous 
system, the second term, gi, accounts for the effect of order parameter variations, and 
the last term, (72? corresponds to the Frank elastic energy of a nematic fluid with spatially 
varying director. This last term relates the free energy functional (171) to the well-known 
Frank elastic energy'^' 

Ff = ^J d\ {iri(Vn)2 + 7^2(11 ■ V X n)2 + 7^3(11 x V x n)^} (9) 
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and allows to identify the three Frank elastic constants splay (Ki), twist (-^'2); and bend 
{K3). We note that Eq. ([7]) predicts Ki = K^, whereas experimentally and in simulations 
(see Tabled]), the parameter is usually much higher than Ki. To improve the theory 
in this respect, one would have to include higher powers of Q in the expansion ([6]). 
Minimizing the free energy ([71) with the boundary conditions S = aX z —>■ —00, 5* = 1 at 
z —>■ 00 and fixed director n = no yields the mean-field structure of a planar NI interface 
at fixed anchoring angle cos 6' = uq^z = (hqN) (N being the interface normal): The order 
parameter adopts a tanh profile, S{z) = SqS{z/w) with 

S{t) = ^(1 + tanh(r)) and the width w = 2^o\/l + a(noN)2. (10) 
The mean-field interfacial tension is 

a = ao^l + a{noNy with (Xo = eo/^o- (H) 

The parameter a thus not only determines the elastic constants, but also the strength 
and the direction of the anchoring at the interface. At a > 0, the interface favors planar 
alignment, and at a < 0, it favors normal alignment. 

To study capillary waves, we must allow the interfacial position to vary. We first consider 
a simplified variant, where the director is still taken to be constant throughout the system, 
n = no, and lies in the {x,y) plane (no,^ = 0). For the order parameter, we make the 
Ansatz ^(r) = SoS[{z — h{x, y))/w], where S and w are given by Eq. flTU]) . We note that 
the local surface normal N is no longer constant, but depends on the gradient of the local 
interface position h{x,y), i.e., N oc {—dxh, —dyh, 1). After inserting this Ansatz into the 
free energy functional ([7]), Fourier transforming {x,y) q, and defining q = q/g, we 
obtain 

F = aoA+iy'AlMq)rg'7(q) with 7(q) :=ao(l + a(qno)2)). (12) 

The capillary wave spectrum can be calculated in complete analogy to the simple fluid 
case,I^^l Eq. ([2]), and one gets 

(|Mq)P) = fc,T/7(q)g^. (13) 

Hence this simple approximation already produces an anisotropic capillary wave spectrum. 
It predicts the proportionality l/{\h{q)\'^) oc suggested by our simulation data. Two 
remarks are in order here. First, the analogy to simple fluids is perfect, if the function 7(q) 
is interpreted as a wave vector dependent interfacial tension. Second, 7(q) only depends 
on the orientation of q, not on its absolute value. As a consequence, the capillary wave 
spectrum is scale invariant, all length scales are equivalent, and the interfacial tension is 
predicted to be anisotropic on all length scales. However, we have argued above that this 
is unphysical - the macroscopic interfacial tension must be isotropic. It turns out that 
the problem is caused by the constant director constraint. To resolve it, we must extend 
the theory such that the director is allowed to follow the interfacial undulations. 
To this end, a second approximation has been adopted: The "local profile" approximation. 
The main assumption here is that the width of the interface is small, compared to the 



6 



relevant length scales of the capillary waves. This is of course questionable, since we have 
just seen that capillary waves in simple systems are scale invariant. However, computer 
simulations of various systems^^^"^^' have shown that the concept of separating "intrinsic 
profiles" and capillary waves often provides a highly satisfactory quantitative description 
of interfacial structures. We separate the free energy ([71) into an interface and a bulk 
contribution, F = Fj + Fb- The bulk contribution, 

Fb = r rfz /rfxc/?/|(3 + 2a)(Vn)2 + (3-«)(n-Vxn)2+ (3 + 2«)(nx Vxn)2|, (14) 

Jh(x,y) J ^ J 

accounts for the elastic energy in the nematic fluid. The remaining interface free energy 
Fi vanishes far from the surface. It is evaluated using the assumption that the local order 
parameter profile has mean-field shape in the direction N perpendicular to the interface, 
and that the director variations are slow, compared to the order parameter variations in 
the vicinity of the interface, such that they can be approximated by a linear behavior in the 
interfacial region. The total free energy is then minimized with respect to the director field 
n(r) for fixed interfacial position h{x,y), and for given bulk director no = limz-,oo'n{r). 
Since globally, we still have planar anchoring at a > 0, the bulk director ng lies in the 
{x, ?/)-plane. 

Details on the calculation can be found in Ref. 16. Here we just sketch the main results: 
The final free energy as a function of h{x,y) can be cast in the same way as Eq. f|T2|) . 
However, the wave vector dependent interfacial tension 7(q) now depends on the absolute 
value of q. Expanded in powers of q, it takes the form 

7(q) ^ ao(l + g C3[(qno)2] + q' C^iquof] + ■■■), (15) 

The leading term is isotropic. The capillary wave spectrum, still given by Eq. f[T^ . remains 
anisotropic, but the anisotropy comes in through the higher order contributions to 7(q). 
Thus the improved theory resolves our problem: At large wavelengths, the interfacial 
tension becomes isotropic and assumes the value 7 = ctq. This is also the surface free 
energy per area that one would obtain in a macroscopic measurement. '^^1 It is worth 
noting that in the direction perpendicular to the director (q ± no), the Cj in Eq. (ITS!) 
vanish, [^^1 such that 7(q) is equal to ctq at a// wavelengths. In that direction, the capillary 
wave spectrum corresponds to that of a simple interface, which is solely determined by 
the macroscopic surface tension 7 (Eq. [2]). 

These results are gratifying. However, the predicted capillary wave spectrum is now in 
apparent disagreement with the simulation data, which seem to point towards a general 
(|/i(q)P) oc 1/g^ dependence. To assess the apparent discrepancy, we must compare the 
theory and the simulation data at a quantitative level. 

Comparison between Theory and Simulation 

At coexistence, the theory has only one dimensionless parameter, a. Since it enters the 
ratios of the elastic constants as well as the interfacial anchoring parameters, we have 
several independent ways of determining its numerical value. In the following, we shall 
estimate a for the ellipsoid system, based on the data collected in Table [1] For the hard 
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Figure 2: Capillary wave spectrum in the direction parallel and perpendicular to the 
director as predicted by the theory for a = 1. The thin solid line corresponds to the 
constant director approximation, Eqs. f|T2l) and f|T3l) . the thick solid line to the local profile 
approximation. In the direction perpendicular to the director, both approximations give 
the same result, which also agrees with the prediction for simple fluid interfaces, Eq. ([2]). 
Left: Linear plot vs. g^. Right: Double logarithmic plot showing the crossover between 
isotropic and anisotropic behavior. 

rod system, the elastic constants and the anchoring parameters are not yet available - for 
particles as elongated as ours, however, the numerical value for a should be comparable. 
According to the data of Ref. 5, the ratio of interfacial tensions for normal and planar 
anchoring is roughly given by TnT"^' IInT" ~ Comparing this with TnT''' HnT" = 
v^l + a, we get a ~ 1 . Two other estimates for a are provided by the numerical values of 
the elastic constants. According to the theory, the ratio of Ki or and K2 is Ki ^^/ K2 = 
(3 + 2a)/ {3 — a). Numerically, Ki and are different, hence we obtain two different 
values for a. From K3/K2 ~ 11.3, we get a ~ 2.3, and from K1/K2 ~ 2.43, we get a ~ 1. 
Comparing the three estimates, we conclude that the choice a = 1. is reasonable. This 
can now be inserted in the theory. 

Fig. [2] shows the corresponding theoretical curves for the capillary wave spectrum. As 
discussed above, the spectrum becomes isotropic in the limit of large wavelengths, q ^ 0. 
However, the crossover occurs in the g-range g^o ~ 10~^, corresponding to the length scale 
27r/g ~ 600^0- Isotropic behavior is expected for wave vectors less than qC,o ~ 3 ■ 10~^, 
i.e., length scales larger than ~ 2000.^o- Taking into account that ^0, is roughly half the 
interfacial width, the theory would thus predict an isotropic capillary wave spectrum on 
the length scale of ~ 7000 particle diameters in the ellipsoid system, or ~ 10000 particle 
diameters in the hard rod system. Hence it is not surprising, that this regime has not 
been observed in the simulations. In the g-range g^o ~ 0.05 — 0.2 (corresponding to 
27r/g ~ 100^0), the capillary waves predicted by the local profile theory follow closely 
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those of the constant director approximation: The director is effectively stiff. The elastic 
penalty on director variations is sufficiently strong that it prevents the director from 
following the interfacial undulations. At large wave vectors g^o ~ 0.4, the curve predicted 
by the local profile approximation drops sharply. This is an artefact, the approximation 
breaks down for such small length scales. ^^^^ We recall that the local profile Ansatz is 
based on the assumption of separate "interfacial" and "capillary wave" length scales, and 
is thus bound to fail as qC,o approaches one. 

Comparing the theoretical capillary wave spectrum. Fig. [21 with the simulation results. 
Fig. [H and disregarding the high g-regime, we find reasonable agreement between theory 
and simulations. The theory predicts that the capillary waves in the direction of the 
director should be suppressed by a factor of the order two, which is roughly reproduced 
by the simulations. It explains why the simulations fail to produce an isotropic fluctuation 
spectrum at large wavelengths, and why the fluctuations in the direction perpendicular to 
the director are rather well described by the simplest capillary wave theory for isotropic 
fluid interfaces, Eq. ([2]), if one inserts the global surface tension for 7. Moreover, it gives 
a reason why the capillary waves in the direction parallel to the director may show a 
behavior even in the anisotropic regime, as has been observed in the hard rod system. 

Discussion and Summary 

In the present paper, we have reviewed and reexamined recent work on interfacial fluctua- 
tions (capillary wave fluctuations) of NI interfaces in situations where the bulk director is 
on average oriented parallel to the interface. In the past, we have studied such interfaces 
by computer simulations of two model liquid crystals, and by an analytical continuum 
theory. Here, we compare the results of the different studies and relate them to each 
other. The interfacial capillary wave fluctuations of NI interfaces are governed by the 
competition of the interfacial tension, the anchoring of the nematic director at the inter- 
face, and the elasticity of the director in the bulk. As a result, they are anisotropic, the 
fluctuations are suppressed in the direction parallel to the director. If the director were 
infinitely stiff, the anisotropy would persist on all length scales. Relaxing the director 
relieves this behavior, and the fluctuations become isotropic for large wavelengths. Hence 
we can identify two regimes: A "flexible director" regime, where the capillary waves are 
isotropic and governed by the macroscopic surface tension, and a "stiff director" regime, 
where they are anisotropic and governed by a "surface tension" 7(q) which depends on 
the orientation and the magnitude of the wave vector. The crossover length scale between 
the two regimes is very large, such that it could not bee observed in the simulations. 
Otherwise, the simulations are in good semiquantitative agreement with the theory. 
We thank Nobuhiko Akino, Michael Allen, Andrew McDonald, Jens Elgeti, and Richard 
Vink for enjoyable and fruitful collaborations that have led to the results presented in this 
paper. This work was funded by the Deutsche Forschungsgemeinschaft (DFG). The simu- 
lations were carried out at the John- von- Neumann computing center at the Forschungszen- 
trum Jiilich. 

Appendix: Fourier transform of h{x, y) 

In this appendix we shall briefly sketch how the data for the interfacial height fluctua- 
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tions h{x,y) have to be Fourier transformed in order to ensure the vahdity of the simple 
expressions ([2]) and (fT3l) . Given a planar interface with the projected area A = L^Ly, 
and a discrete data set h{xn,ym) with points {xn,ym) on a regular grid (n e [l,iV],m e 
[1, M], x„ = L^n/N, ym = Lym/M), then the Fourier transform h{qk, qi) is defined by 

Hqk,qi) = ^Ee^^''='"^''""^M^n,Z/™) (16) 

n,m 

with g-values qk = 2TTk/Lx, qi = 2711/ Ly and k G [1,A^],/ G [1,M]. The inverse Fourier 
transform is h{xn,ym) = Lfc.z g/). We note that h{qk,qi) has the 

unit of a squared length. The convention (fT6!) differs from the symmetrical definition 
^(O'fc, qi) = ^Ijj^j Z]n,m e**-'^'=^"'''^'^'"''/i(a:„, which is often used for discrete Fourier trans- 
forms. When choosing the latter, one has to introduce an additional prefactor A/NM in 
Eqs. (12]) and ([T^. 
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